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Abstract 

The HBT method is used to get information about the sizes, shapes 
and sometimes also about the time evolution of the homogeneity re- 
gions in hadroproduction processes. Homogeneity region K is the 
region, where the hadrons with momentum K are produced. The 
shape and size of homogeneity region K is described by the Wigner 
function W(K, X) evaluated in the interaction representation after all 
the hadrons had been produced. Additional information about the 
evolution in time is contained in the emission function S(K,X). A 
theorem is presented and discussed which specifies which of the pa- 
rameters characterizing the Wigner function can and which cannot be 
measured using the HBT method. In order to obtain the complete 
Wigner function additional assumptions are needed. For instance, it is 
enough to know the distribution of the centers of the homogeneity re- 
gions (X)k- In order to find the emission function further assumptions 
are required. No systematic analysis is available, but some instructive 
examples are discussed. 
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1 Introduction 



The HBT method is a somewhat controversial [T] name for the use of Bose- 
Einstein correlations to get information about the interaction regions, defined 
below, in multiparticle production processes. For a recent review see [2J. In 
spite of its long history since the seminal paper of the Goldhabers Lee and 
Pais [3], the method remains a field of active research. In the present paper 
we discuss what can and what cannot be measured using the HBT method. 
We also describe some open problems. 

The interaction region is interpreted as follows. Consider a high energy 
collision of two heavy nuclei where many hadrons are produced. The tra- 
jectory of each hadron begins somewhere - in the point where the hadron 
got created. This is a classical picture. It can be improved by replacing the 
points of creation by small regions in space. All these points, or regions, 
averaged over many similar collisions, form the interaction region. Thus, the 
interaction region is, in general, not the volume where the fluid (strongly in- 
teracting quark-gluon plasma?) created in the collision evolves. For instance, 
in some models it is a two-dimensional shell. 

If we were not constrained by the laws of Nature, the best description of 
the interaction region would be the time dependent probability distribution 
in phase space, i.e. the probability density for a particle of momentum p 
to be created at space-time point X. Since, however, according to quantum 
mechanics it is not possible to measure simultaneously and precisely the po- 
sition and the momentum of a particle, one must compromise. In ordinary 
quantum mechanics, where the number of particles is conserved, the Wigner 
function is considered to be the best replacement for the phase space prob- 
ability distribution. In hadroproduction the situation is more complicated, 
because the number of hadrons increases from zero to some final number. 
The analogue of the Wigner function applicable in this case is known as 
the emission function. It is implicit already in the work of Shuryak 0], but 
has been first explicitly defined and used by Pratt [5]. For a fairly recent 
discussion see [2J. 
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2 Emission function 



The emission function can be defined^ by the formula (cf. e.g. [2]) 

S(P,X)= Uy Y^n(X + hT R (X-^)e^ Y , (1) 

J R 

where 

X = ^(x 1 +x 2 ), Y = Xl -x 2 . (2) 

T R {Z) is the probability amplitude for producing a hadron at space-time 
point Z. This, of course, depends on the state of the surrounding, denoted 
here by R, of Z . By analogy with thermodynamics, the averaging over R is 
done on the product T^T R and not on the single probability amplitudes T R . 
The integrand depends on two moments of time: X° ± \Y° . Sometimes (see 
e.g. |5J) incoherence in time of the production process is assumed: 

J2Tr(X + ^)T R (X - |) = 5(Y°)T,(X, Y). (3) 

R 

Then there is only one time, occurring on both sides of equation ([1]), and the 
emission function becomes closer to the Wigner function (see Section 4). 

Note that in (CQ) there is neither integration nor differentiation with re- 
spect to the components of p. Thus, momentum appears only as a parameter. 
This is related to the fact, discussed in the following section, that it is not 
possible to measure the whole interaction region. The most one can hope for 
is to measure the homogeneity regions. Homogeneity region K is the region 
where the hadrons with momentum K got created. Thus, the problem is to 
find the profiles of the homogeneity regions^: 

/+oo 
dX°S(K,X). (4) 
-oo 

This limitation of the HBT method was noticed by Bowler [6j. The name 
homogeneity region was introduced by Sinyukow [7]. 



order to make the emission function closer to the Wigner function, we have changed 
the sign of the exponent with respect to the one given in [2] 

2 Here and in the following ~ means: equal up to a known proportionality constant 
irrelevant for our discussion. 



3 



The emission function is related to the single particle density matrix in 
the momentum representation by the formula 



The density matrix does not depend on time, because it refers to a time when 
all the hadrons are already present and propagate freely (To some extent one 
can include final state interactions see e.g. [2]). It is written in the interaction 
picture. 

Formula (jSJ) rises two problems. How to measure p(K, q)? The diagonal 
elements are given by the single particle momentum distribution, but the 
textbook advice to look for distributions of other measurable quantities in 
order to find the out of diagonal elements is inapplicable here, because mo- 
menta are all we know how to measure. A brilliant partial solution to this 
problem [3] is discussed in the following section. Since Wigner's function 
is the Fourier transform of the density matrix, knowing the density matrix 
is enough to find the size and shape of a homogeneity region. The emission 
function, however, contains additional interesting information about the time 
evolution of the hadronization process. 

The second problem is, how to solve for S equation (jSJ) for a given density 
matrix? At first sight it might seem that it is enough to invert the Fourier 
transformation, but in order to do that one would have to know p(K, q) for 
all the four-vectors q and not only for all the three-vectors q as is the case. 
Here not much is known. We discuss the problem in Section 4. 

An important assumption concerning the emission function is the smooth- 
ness assumption, see e.g. [2]. According to this assumption, the dependence 
of S(K, X) on K is so weak that we can replaced K by p\ or p 2 , or replace K , 
by a/ m 2 + K 2 , without changing significantly the results. With this assump- 
tion many objections can be explained away. For instance, two apparently 
very different versions of the HBT method can be shown to be equivalent [2] , 
or the question can be answered: why the arguments of the emission func- 
tion, which are half sums just like the arguments of the Wigner function, 

3 Care is taken to minimize the error. E.g. the product S(K, xi)S(K, X2) is replaced 
by S(pi, X\)S(p2, X2), so that the errors of the two substitutions partly cancel see [2]. 





where 




(6) 
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can be interpreted as particle momentum and position? The smoothness as- 
sumption should hold for pairs of momenta "P\iVi important for the analysis. 
For high energy collisions of heavy ions the relevant momentum differences 
are small and the assumptions seems justified; for e + e~ annihilations or pp 
scattering they are much bigger and the assumption is doubtful [8]. One 
of the outstanding open problems is, why the HBT method is applied with 
comparable success to heavy ion, e + e~ and pp collisions? 

3 Measuring the density matrix 

Much information about the single particle density matrices can be obtained, 
under certain assumptions, by studying the distributions of momenta for sets 
of n = 1,2, .. . identical mesons, for instance n~ mesons |3j. The formulae 
used to relate the n-particle momentum distributions P(p 1 , ■ ■ ■ ,p n ) to the 
single particle density matrix elements are: 

^i(Pi) = N lP ( Pl ; Pl ), 
^(Pi,P 2 ) = P(p 1 )P(p 2 ) + iV 12 |p(p 1 ;p 2 )| 2 , 

PsiVi, P2, Ps) = P(Pi)P(P2)P(P3) + ^(Pi, P 2 ) J P(P 3 ) + (7) 

p (p 2 > p 3 ) p (pi)+ p (Ps> Pi)^(p 2 ) + Ni 23 9fc[p(Pi; p 2 )p(p 2 ; p 3 )p(ps; Pi)] 

and so on, where 3? stand for real part of and N a are normalization con- 
stants, irrelevant for our discussion. The first formula follows from the def- 
inition of the density matrix. The others are derived like the second one, 
which was obtained, for a specific model and in a different notation, in [3]: 
For two totaly uncorrelated particles, the probability distribution for their 
momenta would be proportional to p(pi, p 1 )p(p 2 ; p 2 ); when the two particles 
are identical and have spin zero, symmetrization introduces the correction 
term p(p x ; p 2 )p(p 2 ; Pi) = \p(Pi, P 2 | 2 - 

There is a problem with the consistency of equations (J7J). Suppose that 
exactly two particles are produced. Then F 2 (p 1; p 2 ) integrated over p 2 should 
give Px(Pi)- Actually, the integral of the first term on the right-hand-side 
of the formula for P 2 yields Pi and the second term supplies an unwanted 
correction. The dependence of this correction on p l is usually different than 
that in Pi(p 1 ). Thus, it is not possible to compensate it by a normalizing 
factor. The only way out is to make this correction negligibly small. Since 
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p(Pi> Pi)i s known and fixed, one must assume that p(Pi, p 2 ) decreases rapidly 
with increasing |q|. 

The model described here has a number of other difficulties which, how- 
ever, can be relieved by a suitable handling of the experimental data and/or 
of the comparison between theory and experiment. We will describe them 
very briefly. A much more detailed description with many references can 
be found in the review [2]. The model assumes no correlations except the 
Bose-Einstein correlations. The remedy is to construct a sample which con- 
tains all the correlations except the Bose-Einstein correlations and compare 
it with the full sample which contains all the correlations including the Bose- 
Einstein correlations. The relative momentum q of two charged particles 
with similar momenta is strongly affected by Coulomb, and in some situa- 
tions also by strong, interactions. This is the famous problem of final state 
interactions and is usually handled by introducing a suitable q-dependent 
correction term. Many final hadrons are secondaries originating from decays 
of resonances. The short-lived resonances are no problem, but the long-lived 
ones mimic much bigger interaction regions. Actually, they produce in the 
ratio p^p^p^p ) a narrow peak at small |q| which is below experimental reso- 
lution and is usually corrected for by changing the normalization factor iV 12 . 
Finally, there are purely experimental problems like correcting for momen- 
tum resolution, particle misidentification etc. We assume in the following 
that the consistency condition is satisfied and that the experimental data 
have been fully corrected, so that relations ([7]) can be applied. 

The obvious question is, can one solve equations ([7]) for the single particle 
density matrix p(px;p 2 )? The answer is negative. A simple calculation [9], 
[TU] shows that the predicted momentum distributions (for n=l,2,...) do 
not change when p(p x ; p 2 ) changes as follows: 



where /(p) is an arbitrary real function of a single momentum. This group of 
transformations includes all the modifications of p(p x ; p 2 ) which affect none 
of the momentum distributions. The first ambiguity is not much of a prob- 
lem. It corresponds to the space inversion of the homogeneity region. If the 
homogeneity region is symmetric this has no effect. If it is not symmetric, 
one usually can find a physical argument to chose one of the two possibili- 





and/or 



(8) 
(9) 
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ties. The second ambiguity, however, is much more serious. Let us take an 
example. Choosing 

m = \ £ a ^ ( 10 ) 

i=x,y,z 

where 04 are arbitrary real numbers, one finds that this implies the modifi- 
cation 

S(K,X) -> S(K,X + a(K),X°), (11) 

where the components of a(K) are Oj-fQ. Thus, every homogeneity region, 
except the K = one, gets shifted. Since the possible shifts are a class 
of functions of K, the relative positions of the homogeneity regions can be 
almost arbitrarily changed. This, incidentally, is a proof that at best one can 
hope to measure the sizes and shapes of the individual homogeneity regions. 
Putting all the homogeneity regions on top of each other, one could get a 
lower bound for the size of the overall interaction region, but there is no 
upper bound. 

The previous example is a very special case of a general theorem stating 
what can and what cannot be measured using the HBT method [11]. The 
formulation uses cumulants, so let us quote their definition. For the profile of 
any homogeneity region Pk(x) we can define its characteristic function (e lt ' x ) 
which is, of course, a function of the vector t. Consider now the power series 
expansion 

iog<a->= £ ^^ yyy , («> 

HX- '"V '"Z" 

n x ,n y ,n z J 

where the summation is over three sets of all non-negative integers. The 
coefficients K(n x ,n y ,n z ) are the cumulants of the probability distribution 
Pk(x). The number n x + n y + n z is the order of the cumulant. The cumulant 
is even (odd) when its order is even (odd). 
The theorem consists of three points: 

• Every even cumulant of Pk( x ) can be measured. 

• No odd cumulant of Pk( x ) can be measured. 
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• When the first order cumulants, i.e. (X)k, are given, they fix the 
phase /(p) up to an irrelevant constant and, for (X)k 7^ 0, eliminate 
the freedom of inversion. Thus, the full Pk(x) become measurable. 

Let us consider a few applications of this theorem. The second order 
cumulants are elements of the covariance matrix ((xj — (xi))(xj — {xj)))k 
where Xi = x,y,z. Experimentalist routinely measure the HBT radii and use 
them to calculate the elements of the covariance matrices of the homogeneity 
regions. The theorem shows that this interpretation is not affected by the 
ambiguity Q. In the imaging method (see [12] and references given there) 
one determines from the measured momentum distributions the probability 
distributions J>k( x i — x 2) f° r pairs of points within a homogeneity region. 
This is much less information than given by the profiles of the homogeneity 
regions. For instance, the distributions p K are not sensitive to the averages 
(X) K - However, one finds that the distribution p K depends only on the 
even cumulants of the distribution pk and that it depends on all of them. 
Therefore, it is reliably measurable and it gives much more information than 
the measurements of the HBT radii. 

On the other hand, the centers of the homogeneity regions {X)k are first 
order cumulants and cannot be reliably measured. This is a much shorter 
derivation of the conclusion from example ( flOl) that the relative positions of 
the homogeneity regions are not measurable. 

Of special interest is the third point. The collisions of heavy ions are usu- 
ally described using rather simple-minded tools, like Euler's hydrodynamics. 
It is much easier to believe that such analyses give reasonably the centers 
of the homogeneity regions than that they reproduce all the intricacies of 
their shapes. According to the third point, however, once the distribution of 
the centers is know, everything else can be unambiguously measured without 
further model assumptions. 

4 Emission function and Wigner function 

Much less work has been done on the ambiguities in the solution of equa- 
tion ([5]) for p(K, q) known. Instead, one usually derives from a model, or 
guesses, an emission function and checks whether it yields the momentum 
distributions in agreement with experiment. Since a general analysis has not 
yet been performed, we present here only some partial results. 
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Let us assume that in some reference frame all the hadrons got created 
simultaneously at time t = 0. Then the emission function can be written in 
the form 



S{K,X) = 6(t)S(K,X), (13) 

where the function S(K, X) should be found from equation (j5J). Substituting 
into (j3j), one finds, a unique solution 

S{K,X) ~ S(t)W(K,X), (14) 

where W denotes the Wigner function. Thus, in this case the emission func- 
tion is as good as the Wigner function. A solution of this type always exists, 
assuming that the hadronization process is over at t — 0, but usually it is 
not the realistic solution we are looking for. 

Assuming (in some reference frame) no coherence in time, one can relate 
the emission function to the Wigner function by the formula [13] 



S(K,X)~ W ^ X \ (15) 

W dt is the Wigner function for the particles produced in the time interval 
dt around t = X°. This Wi erner function is normalized to the number of 
particles produced in the time interval dt and not to one as is usual for 
Wigner functions. To mark this difference we introduced the overline. With 
this normalization Wdt tends to zero with dt — > and it has to be divided 
by dt to give a finite result. Note that this solution for S(K, X) does not 
depend on K°. The equation for S(K, X) can by obtained by integrating 
both sides of (j!5p over time. The result is 

W(K,X) ~ Jdt S{K,X). (16) 

Solution (fT5|) is one among the infinity of solutions of (fl6|) and the HBT 
method gives no hint how to find it. 

However, (TIB]) can be used as a sum rule for S(K,X). If assumptions 
are made about the time dependence of S(K, X), it is possible to get one 
parameter of the assumed time distribution from experiment. An example 
is contained in [14]. This paper proposed a once very famous formula to 
measure both the size and the life-time of the interaction region by the HBT 
method. The formula was used for years by many experimental groups. The 
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authors assumed that all the particle sources are produced simultaneously 
at t = 0, that there is no coherence in time and that each source decays 
exponentially with the same life time. Then, in agreement with our analysis, 
they were able to determine this life time of the sources. 

Let us consider a generalization of the previous case. Suppose that 
S(K,X) does not depend on q 2 . This is more general than the previous 
assumption that is does not depend on K°. For instance, it is enough to as- 
sume that, in the spirit of the smoothness assumption, K in S(K, X) can be 
replaced by a/ m 2 + K 2 . Applying the Fourier transformation to both sides 
of equation ([5]) we convert the density matrix into the Wigner function and 
using the identity 



In the K = frame one recovers formula f fT6l) . Actually, the use of this 
frame has a number of advantages (see e.g. [2]) and should be strongly 
recommended. 

5 Conclusions 

It is not possible to measure the size and/or shape of the full interaction 
region using the HBT method. Therefore, efforts concentrate on measure- 
ments of the homogeneity regions. From the classical point of view, the best 
description of the size and shape would be a probability distribution for pro- 
ducing a particle with momentum K at space-time point X. As known from 
quantum mechanics, however, this is not possible. The usual translation into 
quantum physics is to look for the Wigner function W^K, X). When written 
in the interaction picture and after all the hadrons have been produced, this 
Wigner function does not depend on time. Its interpretation as a probabil- 
ity distribution ]?k(X) = VF(K,X) is only an approximation, e.g. because 
Wigner functions can be negative in certain regions. It seems, however, that 
at least for high-energy heavy ion scattering this approximation is reasonable 
(for a recent discussion see [15]). 



Kq = = K^q - K • q, 
q°, we find after trivial integrations 



(17) 



to eliminate 




(18) 
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The theorem presented and discussed in Section 3 states that: The even 
cumulants of pk_(X) can be measured using the HBT method. This, in par- 
ticular, justifies the measurements of the HBT radii and the imaging method. 
The odd cumulants cannot be measured, which explains, in particular, why 
one has to study the separate homogeneity regions instead of the whole in- 
teraction region. The theorem also shows that when the distribution (X)k 
is known, it is possible to measure the full Wigner functions for all the ho- 
mogeneity regions and consequently also the full interaction region. 

In order to learn about the time evolution of the homogeneity regions 
one uses emission functions. Even when the Wigner function is known, the 
evaluation of the emission function requires additional assumptions. When 
these assumptions are too strong, one obtains no more than the input. This is 
illustrated by the example leading to formula (fT4"|) . With a judicious choice of 
assumptions, however, one can obtain additional information. An example is 
the model of Kopylov and Podgoretskii [13] , where the life time of the hadron 
sources can be measured. Of course, such results are no more credible than 
the assumptions used to derive them. 

Nowadays people usually propose full models of hadronization, which 
among other things predict also the results of the HBT measurements. The 
comparison with these experimental measurements is used to test the model 
and sometimes also to fix some of its parameters. The credibility of the 
model, and consequently of the description of the interaction region that it 
offers, is based on a comparison of its prediction with a variety of experimen- 
tal data and not just on the study of the HBT results. The question: how 
much one can learn about the interaction region in a (as nearly as possible) 
model-independent way? remains, however, interesting. 
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